Universal critical behavior in single crystals and films of YBa 2 Cu307_,5 
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We have studied the normal-to-superconducting phase transition in optimally-doped 
YBa2Cu307_i in zero external magnetic field using a variety of different samples and techniques. 
Using DC transport measurements, we find that the dynamical critical exponent z = 1.54±0.14, and 
the static critical exponent v = 0.66 ± O.f for both films (when finite-thickness effects are included 
in the data analysis) and single crystals (where finite-thickness effects are unimportant). We also 
measured thin films at different microwave frequencies and at different powers, which allowed us 
to systematically probe different length scales to avoid finite-thickness effects. DC transport mea- 
surements were also performed on the films used in the microwave experiments to provide a further 
consistency check. These microwave and DC measurements yielded a value of z consistent with the 
other results, z = 1.55 ± 0.15. The neglect of finite-thickness, finite-current, and finite-frequency 
effects may account for the wide ranges of values for v and z previously reported in the literature. 

PACS numbers: 74.25.Fy, 74.25. Dw, 74.72.Bk 



I. INTRODUCTION 

The high critical temperatures, large penetration 
depths, and short coherence lengths of high-temperature 
superconductors make it possible to measure critical fluc- 
tuations in these materials, in contrast to conventional 
superconductors [H, 0. In spite of nearly two decades 
of work, however, there is no experimental consensus on 
the critical exponents of the superconducting phase tran- 
sition in zero magnetic field. By performing both DC and 
microwave measurements on thin films, and DC measure- 
ments on single crystals, and doing careful analysis of the 
data that properly accounts for finite size, current, and 
frequency effects, we are able to provide consistent values 
for the exponents . 

There are two fundamental parameters which charac- 
terize a second-order phase transition such as the su- 
perconducting to normal transition The first is the 
temperature-dependent correlation length, £(T), which 
close to the transition temperature T c varies as 



~ \T/T C — 1| 



(1) 



where v is the static critical exponent. A second param- 
eter is the relaxation time t(T), which close to T c varies 
as 



r~e~ \T/T C -1\~ Z », 
where z is the dynamic critical exponent. 



(2) 
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It is generally accepted that, theoretically, v m 0.67 in 
a superconductor in zero magnetic field, since the phase 
transition belongs to the three dimensional (3D) XY uni- 
versality class pf. The theoretical situation for the dy- 
namical exponent z is less certain. Fisher, Fisher, and 
Huse [lj argue that the number of Cooper pairs is not 
conserved, so that model A dynamics which give 
z = 2, should apply. Other theoretical considerations 
yield z = 1.5 [4, similar to model E dynamics. Lidmar 
Q and Weber [y] present Monte Carlo simulations that 
suggest z w 1.5. 

The exponent v can be determined experimentally 
from a number of static experiments. In zero field, 
measurements of penetration depth, M, B magnetic 
susceptibility, [Tol. [Tl1| specific heat[fjl1l2| and thermal 
expansivity 13] largely agree that the static critical expo- 
nent v ~ 0.67, and indicate that the phase transition in 
zero field belongs to the 3D-XY universality class. (Note, 
however, that there are some measurements which yield 
different results. 0, G3 

In principle DC conductivity measurements, which de- 
pend on both the statics and the dynamics of the order 
parameter near T c , can determine both the static critical 
exponent v and dynamical critical exponent z. The expo- 
nents v and z are expected to be universal, but values ex- 
tracted from conductivity measurements are not consis- 
tent. For example, DC conductivity measurements yield 
a wide range of values for critical exponents: v = 0.63 to 
1.2 and z = 1.25 to 8.3. [H M, E3, El, 0, H HH 

AC measurements can determine both the real and 
imaginary parts of the fluctuation conductivity, provid- 
ing another probe of critical dynamics. [J [H, 0, [23| 
Measurements over a broad frequency range allow one 
to probe the dynamical behavior of the system and di- 
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rectly measure the fluctuation lifetime. [l8T] These exper- 
iments are difficult and seldom done, and the available 
results are inconsistent, with values of z ranging from 2 to 
5.6. Booth et al. investigated the frequency-dependent 
microwave conductivity of YBa?_ Cu^Oi a ( YB C O ) films 
above T c and obtained z = 2.3 to 3.[18| Nakielski et al. 
measured the conductivity of YBCO at low frequency(<2 
GHz) and obtained z « 5.6. [2J| Osborn et al. did a simi- 
lar experiment on Bi^Sr-iC aCuiO%+& and obtained z « 
2.[25j For an optimally doped La2- x Sr x CuOi(LSCO) 
film, Kitano et al. found that their data were consis- 
tent with the 3D-XY model with diffusive dynamics, 
v k 0.67, z « 2 in a certain temperature range 26] . 

Although the critical exponents v and z should be uni- 
versal, we see that there is at present no consensus in the 
literature as to their values for the zero-field transition 
in the cuprate superconductors. 

In this paper we report the results of a variety of com- 
plimentary experiments which yield independent deter- 
minations of the critical exponents. In section II, we 
discuss DC transport measurements on both thin films 
and thick single crystals. We present different ways to 
infer each exponent from measurements, and show how 
finite-thickness effects in the films can confuse interpre- 
tation of the data in the limit of small currents. We also 
show that choice of the proper range of current can avoid 
the finite-thickness effects in films, and show how appli- 
cation of a small magnetic field allows the determination 
of v in both crystals and films. 

In section III, we discuss microwave measurements on 
thin films. Just as with DC measurements, microwave 
measurements require a non-zero current density. How 
the applied microwave current density affects the mea- 
sured response has not been systematically addressed. 
Recently Sullivan et al. argued that a finite thickness 
effects at low current density was the reason for previous 
inconsistent results in DC measurements. The ques- 
tion of whether a finite-thickness effect influences the AC 
measurement and the extracted critical exponents, as in 
DC measurements, inspired us to study the power depen- 
dence of the microwave fluctuation conductivity. We find 
that several length scales play a role in AC conductivity 
measurements, and only after their effects are properly 
accounted for can the underlying critical dynamics be 
understood. As a further check, after completing the mi- 
crowave measurements, we re-patterned the same sam- 
ples and performed DC measurements. 

When finite-size effects are properly accounted for, we 
find that all of our results are consistent. As discussed in 
detail below, we find v is 0.66 ±0.10 and z is 1.55 ±0.15. 



II. DC MEASUREMENTS ON SINGLE 
CRYSTALS AND FILMS 

Our crystals are grown by a flux method using Y2O3 
crucibles to ensure crystal purity [U [28j]. The films 
are prepared by the pulsed-laser deposition technique 
at 850 °C and 150 mbar oxygen pressure on SrTiOs 
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FIG. 1: Schematic plots of electric field E versus current den- 
sity J. (a) E-J plot in log-log scale and (b) dlog(£')/dlog(J) 
vs. J in semi-log scale. 

substrates. [29l [301] Transport measurements were carried 
out using the standard four-probe method. The currents 
were applied along the ab plane for the films and along 
the a-axis for the crystals. All connections to the sam- 
ple are made through double-T low pass filters to reduce 
noise [3l[ . We measure the samples inside a cryostat cov- 
ered with a /i-metal shield so that the residual magnetic 
field is less than 2 x 10~ 7 T. At 96 K, the resistivity of 
the crystal is around 70 fiil-cm, and the resistivity of the 
film is around 90 /ifi-cm. 

In theory, log-log plots of electric field E vs. current 
density J isotherms above T c have positive curvatures 
and display nonlinearities in the high currents, as shown 
schematically in Fig.QJa). 

As T c is approached from above, in the limit of J — > 

a, 

(3) 

Thus, in a log(-E) vs. log(J) plot, isotherms above T c 
exhibit ohmic behavior - a slope of one - at low currents. 
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FIG. 2: (a) The electric field E versus current density J for 
an untwinned YBCO single crystal. The dashed line sepa- 
rates the normal and superconducting phases. The isotherms 
above the dashed line exhibits ohmic behavior at low current 
density. The isotherms below the dashed line display the van- 
ishing linear resistivity. Isotherms are separated by 3 mK. (b) 
Derivative plot of the data in Fig. EJa) . The dashed line sep- 
arates the normal and superconducting phases. The crossing 
of the dashed line with isotherm of 93.839 K is due to the 
noise. Both plots indicate z ~ 1.5. 



The isotherms below T c have negative curvatures and 
display vanishing linear resistance (R — > as J — > 0). At 
T = T c , the critical isotherm is expected to show a power 
law behavior [l[ 



E 



(4) 



which is a line with a slope greater than one on a log(-E) 
vs. log (J) plot. 

In Fig. HJa) we show selected E vs. J curves in a log- 
log plot for an untwinned YBCO single crystal. From 
the figure, all isotherms above the dashed line (trian- 
gles) have positive curvatures and have ohmic response 
in the limit of zero current, and all isotherms below the 
dashed line (squares) have negative curvatures and dis- 
play vanishing linear resistivity. We can verify this cur- 
vature by fitting the E — J curves to a second-order form 
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FIG. 3: Power law fit to Eq. ([8} of the ohmic response of the 
isotherms just above the transition temperature T c . Because 
Rl oc (T/T c - l)" (z_1) , we can use this line to find v. By 
setting T c = 93.837 ± 0.003 K and z = 1.5 ± 0.20, we get 
v = 0.71 ±0.30. 



of log(-E) = ao + a\ log(J) + a2[log(J)] 2 , where we use 
the sign of a 2 to indicate the curvature of each isotherm. 
We find that a 2 is positive above 93.838 K and nega- 
tive below 93.836 K [13]. Thus, from Fig. [1(a) we find 
T c = 93.837 ±0.003 K. The dashed line in Fig.^a) sepa- 
rates the superconducting and normal states of the sam- 
ple and the high-current part can be fit to E ~ ji-22±o.io 
From Eq. ((U (using D = 3), we find 



1.44 ±0.2. 



(5) 



Another way to evaluate the power law behavior is the 
1 " : ~ E ) T vs. J plot [29]. From Eq. g]), at the transition 



V d log J I T 

temperature T c , 



'd\ogE\ 
v d log J / T c 



z + 1 
D-l 



(0) 



thus, the critical isotherm is a horizontal line parallel to 
the J axis in a derivative plot. The critical isotherm 
separates the monotonically increasing isotherms above 
T c and the monotonically decreasing isotherms below T c 
in the schematic Fig. [Tfb). The derivative plot generally 
displays the phase transition more clearly than a basic 
plot of E vs. J. 

In Fig. [2](b), we show a derivative plot of the data 
shown in Fig.[2~](a). The dashed line in Fig.^b) separates 
the normal and superconducting phases. The intercept 
of the dashed line is 1.25 ± 0.10 and is expected to be 
(z + l)/2 from Eq. ©. From the derivative plot, we find 



z = 1.50 ±0.20, 



(7) 



which is nearly identical to the result obtained from From 
Fig. dja), using Eq. (gj). Fig. (b) also qualitatively 
shows the change in sign of a 2 discussed above. 

The exponent v can be found from the low-current 
ohmic behavior Rl above T c by combining Eqs. JTJ) and 
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FIG. 4: (a) T c - T m vs. fi H of an untwinned YBCO single 
crystal. Here, as T c — T m (H) ~ H 1 / 2 " , we can find v from this 
line without assuming a value for z, and find v = 0.68 ± 0.10. 
The inset is the melting line for the crystal up to 7 T. (b) A 
similar plot of T c - T g vs. (j, H for a YBCO thin film (d w 150 
nm). From this curve we find v — 0.63 ± 0.10. The inset is 
the glass transition line for the film up to 6.5 T. 

© Q 

R L oc (T/Tc-iy^. (8) 

The slope of the log(i? L ) vs. log(T/T c - 1) plot in Fig. [3 
combined with Eq. ((8]), determines v. We find 

v = 0.71 ±0.30 (9) 

from Fig. [3] 

We can apply a perpendicular magnetic field and look 
at the transition in finite field. According to Ref. [l[, the 
difference between the critical temperature T c and the 
melting temperature T m i g \ (H) is 

Tc-T^H)-!! 1 / 2 ", (10) 

where v is the zero-field static exponent. Eq. (fTU|) is ex- 
pected to be true for clean crystals, where the transition 



FIG. 5: (a) E — J isotherms for a 150 nm YBCO film in zero 
magnetic field. The spacing between isotherms is 50 mK. (b) 
The derivative plot of some selected isotherms from (a). At 
low current density regime, the phase transition is obscured 
by finite-size effects. The spacing between isotherms is 100 
mK. 

is a first order melting (m) transition, as well as for disor- 
dered films, where the transition is a glass (g) transition. 
We show T c - T m (H) vs. hqH in Fig.HJa) and find 

^ = 0.68 ±0.10 (11) 

from a power law fit. This result is consistent with the 
3D-XY model, and is also consistent with the result ob- 
tained from experiments using Eq. ([5]). 

In Fig. EJa), we show the E — J curves for a 150 nm 
thick YBCO c-axis oriented optimally-doped film. The 
isotherms differ by 0.05 K from 92.075 K to 91.225 K. 
Unlike Fig. Ufa), we cannot find a single straight line in 
Fig. E^a) that separates the isotherms into two groups 
which are either concave or convex exclusively. To help 
find the true critical isotherm, we again use the deriva- 
tive plot , where the critical isotherm ideally will cor- 
respond to a horizontal straight line, as in Figs. [T{b) 
and^b). However, in Fig. [5jb) , there is no horizontal 
isotherm, and there are isotherms monotonically decreas- 
ing above 2 x 10 7 A/m 2 , and also monotonically increas- 
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ing below 2 x 10 7 A/m 2 . So, if the experimental setup 
were to allow us to measure even smaller voltages, we 
would expect all of these isotherms would bend down to- 
ward 1 (ohmic behavior) in the derivative plot at smaller 
current densities. 

The cause of this behavior is most likely finite-size ef- 
fects [l9|, [33|. Below the transition temperature, thermal 
fluctuations take the form of vortex loops. 34[ As dis- 
cussed in the Appendix, vortex loops with length scales 
of order 



k B T 

27T$ J 



(12) 



are probed by current density J 0, [33[ . The loops with 
length scale smaller than Lj will shrink and cause no 
dissipation. At high current density, such that Lj is 
less than the thickness of the sample d, the vortex loops 
probed in the experiment are still 3D-like. However, at 
low current density, such that Lj > d, the size of the 
vortex loops probed will be limited by the thickness of 
the sample and vortex anti-vortex pairs will be probed. 
This will lead to a non-diverging energy barrier causing 
ohmic behavior even below the bulk transition tempera- 
ture. When Lj is equal to 150 nm, the thickness of the 
film used to produce the data in Fig. [5j the crossover cur- 
rent density is on the order of 5 x 10 6 A/m 2 , which is close 
to where the isotherms bend towards ohmic behavior in 
Fig.Hb). 

Because of the finite-size effects, the conventional 
method picks an incorrect critical isotherm and expo- 
nent, contributing to the inconsistent results from pre- 
vious transport experiments on high-T c films. However, 
high-current data are not affected by finite-size effects, 
and we can extract the T c and z from the high-current 
regime 19]. If we only look at the high-current regime in 
Fig. [5jb) , it looks very similar to the schematic derivative 
plot (Fig. [T|b)) and the actual derivative plot of crystal 
data (Fig. Mb)). The dashed line in Fig. [jjb), which 
coincides in the high current regime with the isotherm 
of 91.825 K, separates the two phases of the film. The 
transition temperature determined from the high-current 
regime is T c — 91.825 ± 0.025 K, and the intercept of the 
dashed line is 1.27 ±0.07. According to Eq. from the 
high-current data, 



z = 1.54 ± 0.14 



(13) 



which agrees with the result from the crystal data. In ad- 
dition, in our other c-axis oriented YBCO films with the 
thickness d ranging from 100 nm to 300 nm, we get consis- 
tent values of z ranging from 1.43 to 1.6 [35f . In passing, 
we note that one should be cautious about the adverse ef- 
fect of joule heating when the measurement is to be done 
in the high-current regime to avoid the finite-size effect. 
In this regard, while the result reported in Ref. [36] is 
interesting in that it was the first I-V measurements of 
high-T c nano-strips, the extracted critical exponents were 
likely to be inaccurate because of the difficulty of avoid- 
ing Joule heating in nano-strips at high currents. In con- 
trast, we have tested heating in our samples by using the 



low- frequency technique of Koch et aL,[37[ and we have 
found that heating does not affect the DC data in sam- 
ples similar to those measured in this paper at current 
densities less than » 10 9 A/m 2 . [H] 

We are not aware of any way to remove finite-size ef- 
fects that will allow us to use Eq. ([5]) to determine v in 
films. Instead, we use Eq. (jTUJ) which is also applicable to 
the vortex-glass transition. By applying a magnetic field 

we introduce a magnetic length scale I b oc y^f 5 which is 

smaller than the film thickness for B > 0.1 T, effectively 
removing the finite-size limitation in the film. We show 
T c - T g (H) vs. H in Fig. gjb) and find 



0.63 ±0.10. 



(14) 



It is important to note that there is more disorder in 
the film than in the crystal. Besides the point-like oxygen 
vacancy disorder as in the untwinned crystals, there are 
other kinds of disorder existing in the film such as twin 
boundaries and lattice mismatch caused by the substrate. 
However, the similar values of z and v for the untwinned 
crystal and the film argue that the universality of the 
phase transition of high-T" c materials is not affected by 
disorder. 

Hence, taking into account results of DC transport 
measurements on both thin films and thick single crys- 
tals, we obtained the critical exponents 



z = 1.54 ±0.14 
v = 0.66 ±0.10. 



(15) 
(16) 



III. 



MICROWAVE MEASUREMENTS ON THIN 
FILMS 



The samples we used for microwave measurements are 
YBCO films (d =100 nm to 300 nm thickness) deposited 
via pulsed laser deposition on NdGaO^ and SrTiO^ sub- 
strates. AC susceptibility showed T c of the films around 
90 K with transition widths about AT C = 0.2 K. The 
resistivity of the films is about 120/iil-cm at 2 K above 
T c . Using a Corbino reflection technique, we measured 
the complex resistivity p = p\ + ip2 of the samples over a 
wide frequency range. The measured complex resistivity 
is converted to conductivity and the mean field contri- 
bution, as determined from the dc resistivity measured 
from room temperature down to the lowest temperature 
in the same experiment, is removed. [l8l l39l [4lll The pro- 
cess is similar to the method described in [18} to obtain 
the fluctuation conductivity at i . (iij 



A. Frequency Dependent Fluctuation Conductivity 
and Power Dependence 

According to Fisher-Fisher-Huse (FFH), in zero mag- 
netic field when the current density is small the complex 
AC fluctuation conductivity should scale as [l| 



a f l(T,uj) w £ 



z+2-D 



S±{ut). 



(17) 
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In Eq. (fT7|) , £ is the correlation length and r is the fluc- 
tuation lifetime. The function S± is a universal scaling 
function above (below) T c , which should be the same for 
all members of a given universality class. As tempera- 
ture approaches T c , both £ and t will diverge according 
to Eqs. (]} and ©. 

The scaling functions behave as S+(y) — > real con- 
stant and S-(y) — ► l/(—iy) for y — * 0, reflecting the 
low frequency behavior above and below T c respectively. 
As y — > oo, representing T — + T c , S+(y) » S-(y) ~ 
gy(J3-2)/(z-i) wnere g i s a complex constant and D is the 
dimensionality of the system. [l|, [22j 

The complex fluctuation conductivity can be written 
as <Tfi = \afi\e 1 ^" , so both the magnitude and phase are 
predicted to scale 

M « r +2 - D is±K z )i, (is) 

<k = $±K Z ) (19) 

where $± is the phase the of the scaling function S±. At 
T c , one expects [22j 

M~c^+ 2 -^, (20) 

and 

^ = ^(z + 2-D)/z. (21) 

Fig. [6] sketches the expected Fisher-Fisher-Huse AC scal- 
ing behavior of the ma gnit ude and phase of fluctuation 
conductivity near Tf.fll |22| 

Fig. [7| shows the measured complex fluctuation con- 
ductivity vs. frequency for various temperatures at two 
different microwave powers. These data display signifi- 
cant and systematic deviations from the expected FFH 
scaling sketched in Fig. [6j At high frequency, both the 
magnitude and phase of the fluctuation conductivity look 
similar to FFH theory. However, as frequency decreases 
the measured magnitude of the fluctuation conductivity 
below T c saturates, instead of bending up. All of the 
phase isotherms below T c tend toward zero, indicating 
ohmic response, instead of approaching tt/2 at low fre- 
quency. These deviations are qualitatively similar to the 
low current-density deviations of E vs. J in DC measure- 
ments seen in Fig. [5] 19] 

Fig. [7] also shows that the applied microwave power af- 
fects the measured fluctuation conductivity, particularly 
at low frequencies. As frequency decreases, the higher 
applied microwave power decreases the magnitude of the 
fluctuation conductivity and depresses the phase. These 
phenomena cannot be explained by the AC scaling equa- 
tion, Eq. (|17[) . and we need to look at the full version of 
the FFH dynamic scaling function, which can be written 
in the following form with assumed dimensionality D = 3 

j=e- Z X±(JfM z ,Hf,-)- (22) 
where E is the electric field. 




FIG. 6: Schematic plots of (a) Magnitude \ofi\ vs. u in 
log-log scale and (b) phase </v vs. lo in semi-log scale at 
various temperatures around T c , based on Fisher-Fisher-Huse 
AC scaling. [J 

Since the critical point is located in the limit of zero 
magnetic field H, current density J and frequency u, in- 
creased applied current should drive the system further 
away from the transition and thus into the ohmic regime. 
In our measurement, the magnetic field term H£ 2 can 
be ignored. The two remaining terms are J£ 2 and w£ z . 
Qualitatively, at low frequency, lu^ z is small so that the 
applied power term, J£ 2 , has more effect on the fluctua- 
tion conductivity. 

To illustrate the effect of different powers, vs. mi- 
crowave power at different frequencies is plotted in Fig. 
[8l The power dependence of \<Jfi\ clearly varies with fre- 
quency. At low frequencies (60 MHz, 80 MHz and 100 
MHz), \<Jfi\ vs. incident power increases first as power 
increases (Region I) and then saturates (Region II). At 
very high power, \<Jfi\ decreases again (Region III). At 
high frequencies (> 0.5 GHz) the fluctuation conductiv- 
ity is almost power independent. 

The important features in Fig. [8] are that large ap- 
plied power affects the fluctuation conductivity, and that 
even small power depresses the fluctuation conductivity 
at low frequency. While the high-frequency and high- 
power data in Fig. [8] are consistent with Eq. (|22|) . and 
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FIG. 7: (Color online) (a) Magnitude |cr/;| and (b) phase 4> a 
vs. frequency at various temperatures for a typical YBCO 
film(xuhl39). The black lines were measured with -22dBm 
power while the red lines were measured with -2dBm at the 
same temperature. (For clarity, only every other isotherm is 
shown.) 
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FIG. 8: (Color online) \ofi\ vs. incident microwave power at 
different frequencies. ( T=89.140 K, sample xuhl39 below 



thus can be explained by FFH scaling theory [l|, [22|, l23j |. 
the low-power low- frequency behavior is not consistent. 

The similarity between this low power and low fre- 
quency deviation and the low current density deviation 
in DC conductivity measurement suggests the presence 
of a " probed length scale" for a finite frequency. As dis- 
cussed in connection with Eq. (|12p and in the appendix, 
when a current with density J is applied, vortex loops 
(with large r) will "blow out" to infinite size (producing 
dissipation). Vortex loops with small r shrink and anni- 
hilate (with no dissipation). A current density induced 
length scale Lj, given in Eq. (| 1 2[) separates vortex loops 
into two categories, depending on their ultimate fate. 

The shrinking of a loop takes time. This time de- 
pends on the size of the loop, thus relating the size of 
a vortex loop to a time scale. In AC measurements, 
small frequency means that large length scales are probed 
and vice versa. By generalizing the order parameter 
relaxation time scale in time-dependent Ginzburg Lan- 
dau theory 40] one can construct a frequency-dependent 
length scale 



<^)'/-«o), 



(23) 



where c is a constant of order 1 and £(0)/£(T) = \T/T C — 
ll". 

In AC conductivity measurements, the probed length 
scale should be determined by both frequency and cur- 
rent density. Since the smaller length scale dominates the 
measured fluctuation conductivity, we propose a plausi- 
ble expression for the probed length scale for AC mea- 
surement Lac, 



1 



1 



1 



Lac Lj L u 



(24) 



This formula has the correct limits as J — > or u> — > 0, 
which corresponds to frequency dependent o>£ 2 scaling 
or current density dependent J£ 2 scaling, respectively. 
It is also qualitatively consistent with the two-term FFH 
scaling without the magnetic term Ht; 2 of Eq. (|22|) in the 
crossover range. Finite size effects come into play when 
Lac approaches the thickness of the film. 

Fig. [5] summarizes the length scales in an AC measure- 
ment in terms of experimental quantities. [43j In this 
figure, we use £(0) = 5A, c = 1 and z = 1.5. The dotted 
line in the figure gives the boundary Lj = L u . To the 
right and below the dotted line, when L w <C Lj, the fre- 
quency induced length scale dominates, and one observes 
mainly frequency dependent scaling of the fluctuation 
conductivity. Above the dotted line, when L w 3> Lj, 
current-induced nonlinear effects will dominate the be- 
havior. This explains the features shown in Fig. [7] and 
Fig. [51 where the current density has less effect on the 
fluctuation conductivity at high frequency and a larger 
effect at low frequency. 

At low frequency and small current density, Lac may 
approach the thickness of the sample (d) or some other 
length scale that interrupts the fluctuation vortex loops. 
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FIG. 9: (Color online) Summary of length scales and finite 
size effects in Corbino AC measurements of fluctuation con- 
ductivity of YBCO films near T c . The dotted line in the fig- 
ure gives the boundary Lj = L u . At low frequency and small 
current density, the probed length scale Lac approaches the 
thickness of the sample. 



Hence deviations from the simple scaling theory are ex- 
pected when Lac > d. 

In our AC measurements, we want to keep to the limit 
Lac < d to avoid finite-thickness effect. Hence we choose 
to stay at low J but high to. In this region we can find 
the true critical behavior without getting into any finite- 
size effect or crossover difficulties. Our previous analysis 
strayed out of this region and this may account for the 
larger values of z reported before [l8[ and elsewhere in 
the literature. 



B. Improved Data Analysis Method 

In this paper, with very small applied microwave 
power, -46dBm (corresponding to J < 2.2 x 10 5 A/m 2 ), 
and high frequency data, we investigated the frequency 
dependent fluctuation conductivity around T c . Conven- 
tionally, examining experimental data with the scaling 
formulas one can search for the temperature at which 
the conductivity magnitude best fits to a power law and 
has a constant value of <f> a , to determine T c and the dy- 
namic critical exponent z. In this analysis process, the 
determination of T c is crucial because it directly affects 
the value of z. Hence we improved the temperature sta- 
bility and conductivity calibration techniques in the ex- 
periment, enabling the measurement of high quality data 
at small temperature intervals (50 mK). 

Using this data, we improved the conventional data 
analysis method [18| to determine T c . One expects a 
power-law behavior of \<Jfi\ on frequency at T c , with a 
change of curvature on either side (a convex function be- 
low T c and a concave function above T c ), as sketched in 
Fig. [|3a) . One also expects a plateau in the conductivity 
phase vs. frequency at T c , with a change in the sign of 



the slope on either side, as sketched Fig. EJb). 

Unlike the DC I-V curve where Strachan et al. used an 
opposite concavity criterion to determine T c in a dL/dV 
plot,[2i| it is hard to take the frequency derivative of 
|<7/;(o;)| because of noise. An alternative approach is to 
do a quadratic fit to the data on a log-log plot. Below T c , 
the curve bends up with a positive coefficient of [log(uj)] 2 
and above T c , the curve bends down with a negative co- 
efficient of \log(uj)] 2 . Hence we did a quadratic fit and 
found that the coefficient of the [log(ui)] 2 term changes 
sign between temperatures 89.192K and 89.245K, brack- 
eting T c . 

The scaling theory also predicts a constant phase angle 
4>o-(u>) at T c . 4> a {uj) vs. logw is known to be a decreasing 
function below T c and an increasing function above T c . 
A linear fit of (f> a (u) vs. log/ also has been done and 
the result shows its has negative slope at 89.192K and 
positive slope at 89.245K, which is consistent with the 
quadratic fit result of log\<7fi(uj)\ vs. logo;. The next 
step is to do a linear fit for log |ct//(w)| verses log(uj) to 
get the slope of log \afi(uj)\, and take the average of the 
4>a(w) at T c to obtain the value of z. From this method, 
we get the critical temperature T c = 89.22 ± 0.05i4T and 
the critical exponent z = 1.62 ± 0.20. 

In addition, we developed a new method to determine 
T c from the data. Consider the Wickham and Dorsey 
scaling function above T c 22\ 



S+(y) 



2z 2 [l 



D-2- 



Hy-{\-iy)( D - 2+z y z ] 



y 2 {D-2-z)(D-2) 



(25) 



where y = lot oc w£ 2 . We find at small y, correspond- 
ing to temperatures far above T c , the function S+(y) is 
essentially independent of dimensionality D and z be- 
cause the fluctuation contribution is small. According 
to Eq. (fTT)) . one can write <?fi(T, to) = ao(T)S(Lu/LUo) 
where <Jq{T) and loq{T) are characteristic conductiv- 
ity and frequency scales, respectively. Both the phase 
</> CT (= tan^ 1 ^!' / 'a{ 1 ]) of <jfi and the magnitude cr/zl/oo 
can be treated as scaled quantities with two temperature- 
dependent scaling parameters coo(T) and ao(T). This is a 
new data collapse method, pioneered by Kitano et aZ.(26j 
They pointed out that the advantage of this new collapse 
method is the independence of the two scaling parame- 
ters Wo (T) and cr (T) ■ In this new data analysis method, 
the parameters u>o(T) and o~q(T) are chosen at each tem- 
perature to collapse <pa(T) vs. lo/loq and \afi\/ao(T) vs. 
lu/luq to smooth and continuous curves, without a priori 
determination of T c or critical exponents. 

First u)q(T) is determined through a collapse plot of 
4> a vs. io/loq{T) from high temperature to low temper- 
ature (see Fig. fTUTa)). Using the feature that S + (y) is 
not sensitive to dimensionality D and z far above T c , 
the appropriate loq(T) for isotherms far above T c is cho- 
sen to make the measured <j) a {uj / ljq{T)) overlap with 
the theoretical prediction from the known scaling func- 
tion 4>(S+(y)). Then at temperatures closer to T c where 
S+(y) starts to depend on D and z, loq(T) for each tem- 
perature is chosen to connect smoothly to the existing 
curve of (j)<j(uj/LUo(T)) and to make all the temperature 
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FIG. 10: Scaling of phase and magnitude of fluctuation con- 
ductivity for sample xuhl39 to determine u>o(T) and <7o(T). 
(a) cj) a vs. u/lo (T); (b) \cr fl \/a (T) vs. lu/lo {T). Solid lines 
are the theoretical calculation of \S+(y)) \ and 2<j)s + (y)/ n from 
Eq. (|25[1 for different values of z, assuming D — 3. These 
lines can be used to determine U}q(T) and <Jo(T) far above 
T c . In this figure, only the uio{T) and co(T) of the isotherm 
T = 89.763 K are shown to make the measured tj>a-{u/uo(T)) 
and \cTfi\/ao(T) vs. ui/u}q{T) overlap with the theoretical pre- 
diction. For all the other isotherms, Wq(T) and t?o{T) for each 
temperature are chosen to connect smoothly to the existing 
curve of /u>o{T)) and |<t/; \/ao{T) respectively, to make 
all the temperature curves collapse into one smooth and con- 
tinuous curve. 



curves collapse into one smooth and continuous curve. 
This process continues to lower and lower temperature 
until a temperature is reached where <j) a (uj /lvq(T)) can 
not be connected smoothly to the existing curve. In this 
way, loq(T) for temperature points above T c can be de- 
termined. 

To scale the conductivity magnitude, we start with 
the determined loq(T) for each temperature, then plot 
\<jfi\/ao(T) vs. lu/loo(T), where <jq(T), similarly to 
u>o{T), is determined for each temperature to make 
a smooth and continuous curve of \afi\/aa(T) vs. 
w/w (T).(seeFig.[HH>)) 

Using the power-law assumption for loq(T) and cto(T), 
T c can be determined. Fig. 1111 shows u>q(T) vs. t and 
cto (T) vs. t for different assumed values of T c . The cor- 
rect T c can be determined from the line showing a pure 
power-law. Fig. [TTT a) shows that the blue line which 
corresponds to an assumed T c = 89.25K is straightest. 
Fig. fTTTb) also shows that the blue line is straightest. 
From these two figures, T c is consistently determined to 
be T c = 89.25 ±0.05/1'. This result is also consistent with 
the T c determined by the improved conventional method. 

With the value of T c determined here, we can do a 
linear fit for log verses logiu) to get the slope 

of log |<7/z(a;)|, and take the average of the ^cr(w) at T c 
to obtain the value of z. Through this procedure, we 
obtained the critical exponent z — 1.55 ± 0.20. 

In the procedure outlined above, we take advantage of 
the broad microwave frequency range of the experiment, 
which includes frequencies of order 1/r. High quality 



FIG. 11: lu (T) vs. [2^1 and a Q (T) vs. for different 

assumed T c for sample xuhl39 and temperature from 89.297K 
to 89.763K. The errors of wo(T') and <Jo(T) are about the size 
of the points. 

data at small temperature intervals are essential for the 
implementation of this method. Another advantage of 
this method is that many isotherms near T c contribute to 
defining the scaling curve, not just the one closest to T c . 
The new method has the advantage of more precisely de- 
termining T c . So according to the two methods the crit- 
ical temperature and exponent for sample xuhl39 were 
determined to be T c = 89.25 ± 0.05^, z = 1.55 ±0.20.0 
The dynamic critical exponent should be sample in- 
dependent. To check the results, we not only repeated 
measurements on the same sample, but also repeated the 
experiment on different samples. Films of different thick- 
ness (d =100 nm to 300 nm) were examined, and z was 
found to be independent of the thickness, keeping in mind 
the constraints of Fig.[5J Experiments on 6 samples have 
been done giving 

z= 1.55 ±0.15. (26) 

C. AC and DC Experiments on the Same Sample 

We also performed DC current-voltage character- 
istic measurements on the same samples. [4lj Typi- 
cal results are shown in Fig. 1121 (with no background 
subtraction 30] ) According to the negative curvature 
criterion[2i|, we determined the critical temperature to 
be 91.220±0.04K and the critical exponents = 1.75±0.2 
from the derivative plot in Fig. [T2Tb). In Fig. [TH all 
the isotherms tend towards ohmic behavior at low cur- 
rent density, brought about by Lj > d finite-size effects, 
as discussed for the data shown previously in Fig. [51 [l9| 
From this data it is clear that when the current density 
is smaller than 1 x 10 6 A/m 2 , the sample will have only 
ohmic response around T c . The -46dBm applied power in 
the AC measurement corresponds to a maximum current 
density of 2.2 x 10 5 A/m 2 (< 1 x 10 6 A/m 2 ). This means 
that for -46dBm incident power Lj > d, verifying a fea- 
ture of Fig. [51 and suggesting that one-parameter scaling 
should work when L u < Lj, d. Hence it is appropriate 
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FIG. 12: (Color online) DC current-voltage characteristics 
measurement, performed after the AC experiment on xuhl39 
in zero magnetic field, (a) E-J isotherms(50 mK apart), (b) 
d log E/d log J vs. J derivative plot(40 mK apart). 

to determine T c and critical exponents with AC data at 
-46dBm applied power. 

The difference of T c between DC and AC measure- 
ments is due to the different thermometer positions and 
temperature control techniques of the two experimental 
systems. The resistance vs. temperature plots from the 
AC and DC experiment have a temperature offset about 
2.0 K, which is the difference of the determined T c from 
these two methods. 

In DC measurements, disorder and heating lead one to 
systematically choose a lower temperature isotherm as 
T c , resulting in an enhanced value of z.[3(| SH We find 
that films with lower normal-state resistivity and smaller 
AT C have smaller values of z.[3(3| The films used in AC 
conductivity measurements were grown on NdGaOs sub- 
strates, and these films have systematically higher resis- 
tivity and larger AT C than films on SrTiO^ substrates. 
In addition, performing DC measurements on the same 
film after AC measurements involves more processing 
steps than a DC measurement alone, and may result in 
additional disorder in the sample. This correspondingly 
gives larger values of z (Fig. [12]) . We carefully repeated 



the DC measurements alone on different YBCO films 
grown on different substrates (SrTiOz and NdGaOz) 
and found that the sample quality does affect the ob- 
tained value of z.[4l[ However, for films with high T c , 
sharp transition and small resistivity, the obtained value 
of z » 1.50, which is consistent with the AC result 
z = 1.55 ± 0.15. In addition, DC measurements car- 
ried out the same way on high-quality crystals shown in 
the previous section also gave z ~ 1.5. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we performed of a variety of complimen- 
tary experiments to determine the critical exponents in 
optimally-doped YBa2Cu307-«5. 

The DC transport measurements on both thin films 
and thick single crystals show how finite-thickness ef- 
fects in the films can confuse interpretation of the data 
in the limit of small currents. Only with the choice of 
the proper range of current, can one avoid the finite- 
thickness effects in films and obtain the correct expo- 
nent z, consistent with the value obtained from thick 
single crystals measurements (where finite-thickness ef- 
fects are unimportant). We also show how application of 
a small magnetic field allows the determination of v in 
both crystals and films. Using DC transport transport 
measurements, we find that the dynamical critical expo- 
nent z = 1.54 ± 0.14, and the static critical exponent, 
v = 0.66 ± 0.10 for both films and single crystals. 

Microwave measurements on thin films at different fre- 
quencies and at different powers have also been per- 
formed, which allow us to systematically probe different 
length scales in the sample. After developing a com- 
prehensive understanding of length scales in microwave 
measurements, we choose to stay at low J but high u) 
to find the true critical behavior without getting into 
any finite-size or crossover effects. DC transport mea- 
surements were also performed on the films used in 
the microwave experiments to provide a further consis- 
tency check. These microwave and DC measurements 
yielded a value of z consistent with the other results, 
z = 1.55 ±0.15. 

To conclude, using two different measurement meth- 
ods, we studied the dynamic fluctuation effects of 
Y Ba2CuzOi-s single crystals and thin films around T c . 
The results of both AC and DC measurements agree with 
the XY value for v w 0.67 and with model-E dynam- 
ics value for z = 1.55 ± 0.15. [|| The neglect of finite- 
thickness, finite-power, and finite-frequency effects may 
account for the wide ranges of values for v and z previ- 
ously reported in the literature. 
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Appendix: Currents and Length Scales in 
Superconductors 

In this appendix we consider a number of length 
scales in current-carrying superconductors to provide a 
clearer physical meaning for the length scale Lj of Eq. 
(fT2")) . We first consider a simple model for fluctuations 
in superconductors, where we assume that the only 
fluctuations are circular vortex loops (or vortex "smoke 
rings") of radius r.[34[ The energy of such a loop can be 
written as 



Ui oop = 2irre(r) (27) 

where e(r) is the energy per unit length of the vortex 
loop. For a straight vortex, 

*-->-£(*)'*G)-£(*H?) 

(28) 

where Kq is a modified Bessel function of the second 
kind and the approximate form holds in the limit of high 
k = X/£ . [45l] As a first approximation, we will assume the 
energy per unit length is constant, given by Eq. (|28p . 

In an infinite superconductor with no applied current, 
vortex loops of different sizes occur with different proba- 
bilities as thermal fluctuations. The probability of find- 
ing a loop of size r in a range dr is given by 



2. t (r) 

P ^ dr = ( 29 ) 

J, 'e k B T 1 dr 

where interactions between the loops are neglected for 
simplicity. 

We wish to find the size of a typical vortex loop, r met j. 
One way to do this is to find the fraction of loops / with 
a radius greater than the median radius, or r > r me d- 
This fraction will be / = i, given by, 



r°° e k B T r ( i r 1 
f°°e-^ r dr ~ 2 ' 



(30) 



If the energy per unit length of the loop is given by Eq. 
|), then Eq. flUD leads to 

r med = ^+^—ln2. (31) 

Z7T£ 



If the second term on the right hand side of Eq. ([31 
dominates, this gives 



W«^ln2^ £ «-MLi n2 . (32) 

Eq. (|3"2"]) states that, within a factor of In 2, the total 
energy of a vortex loop of size r mec i is equal to fc^T, 
which is a plausible result. 
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To check whether the second term on the right side of 
Eq. (f3"Tj) is the dominant one, we combine Eqs. (f2"5| and 
(ED. This leads to 



In 2 



<i> 2 



2tt In/ 



In 2 k 2 
2tt In k 



(33) 

where At is defined in Eq. (1.1) of Fisher, Fisher, and 
Huse [lj (in cgs units with the Boltzmann constant ks 
defined to be 1). The second terms in Eqs. (|31jl and |33|) 
dominate in the critical regime because £ diverges while 
At is fixed. 

For simplicity, we drop the In 2 in Eq. (f31~j) . and use 
the physically plausible result 



Vrned 



k B T 
2tt£ ' 



(34) 



Next consider that a current per unit area J is applied 
in a direction perpendicular to the plane of the loop. The 
total Lorentz force acting outward on the loop is 



F„ 



27rrJ$o- 



(35) 



The energy defined in Eq. (|27|l gives rise to an inward 
force that the loop exerts on itself, — 2-7T£. Summing the 
forces and finding the point where the net force is equal 
to zero leads to a critical loop size 



(36) 



where, for simplicity, e(r) is again assumed to be inde- 
pendent of r. Note that Eq. (|36|) is not the equation for 
the current-dependent length scale Lj. 

Physically, if a vortex loop has r > r 0l the external 
current "blows out" the loop to infinite size; this process 
leads to dissipation. If r < r a , the vortex loop shrinks 
and annihilates. 

One can interpret Eq. (|36[) in a different, but equiv- 
alent, way. The presence of a current density J signifi- 
cantly alters the population of vortex loops with r > r a , 
and has less effect on the vortex loops with r < r D . In 
this sense, a current J probes the physics on length scales 
of order r Q and larger. This is the type of language that 
is sometimes used to describe Lj. 



We next discuss the physical significance of comparing 
the lengths r Q and r med , Eqs. (J34J) and ((36|). If r med < r Q 
, the current is probing a length scale where there are very 
few vortex loops. The current thus acts as a very small 
perturbation on the system. If r mec [ 3> r a , the current 
is probing a very short length scale, and a large portion 
of the intrinsic vortex population is being disrupted by 
the current. The point where r me d = r a thus marks a 
crossover in the behavior from current acting as a small 
perturbation to current acting as a large perturbation. 

How does a non-infinite film of thickness d affect the 
physics? It is plausible to say r me d <C rf is the three- 
dimensional limit, while r me d 3> d is the two-dimensional 
limit, since in the second case most of the vortex loops are 
interrupted by the film thickness, while in the first case 
they are not. This is true as far as it goes, but misses the 
key point that an applied current probes physics at the 
scale of r D and larger. Thus, even in the limit r me d S> d, 
if r Q is small enough, current will probe physics on length 
scales smaller than d, and thus the measurement will not 
be affected by the finite thickness of the film. 

In order for the thickness of the film to have a measur- 
able effect, the current should probe a significant fraction 
of the loop population and should also probe lengths on 
the scale of the film thickness. For this to be true, we 
require 



= r med = Lj 
Combining Eqs. d34j), ([36]), and (J37j) gives 

f k B T \ ^ 



Lj = 



V27r$oJ 



(37) 



(38) 



This argument leading to Eq. (j3"5|) motivates a physical 
description for Lj: For any J there is a length scale Lj, 
given by Eq. (f38|) , such that roughly half the equilibrium 
(zero current) vortex loop population is strongly affected 
by J. This is the length that one should compare to the 
film thickness for seeing whether or not measurements are 
in the two or three dimensional limit. The requirements 
are that there be a significant fraction of the loops that 
would exceed the film thickness, and, in addition, that the 
current is probing the same length scale. 
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